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QJ ' In realistic technicolor models containing many fermions, the electroweak baryogenesis 

i-i^ offers a natural scenario for generating baryon number asymmetry. One of the key ingre- 

dients is the occurrence of the first order phase transition at finite temperature. As a first 
step toward the exploration of this possibility on the lattice, we develop an agile method 
to identify the critical mass for a given Nf , separating the first order and the crossover 
transition. We explain the outline of our method and demonstrate it by determining the 
f0^ , critical mass of A'^j'-flavors in the presence of light two-flavors. It is found that the critical 

^^»^ ' mass becomes larger with Nf. 
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1. Introduction 

Technicolor (TC^Phas been known as one of the most natural candidates for physics 
beyond the standard model. The classical TC models have been already excluded 
for many reasons, while those with many fermion flavors in the fundamental rep- 
^ . resentation are expected to escape from various experimental constraints due to 

special properties called "walking" dynamics.l^ Whether many-flavor TC models re- 
ally show expected walking behaviors is now actively and rigorously investigated on 
the latticcl^l 

In this article, we focus on another aspect of TC models, the nature of the 
thermal phase transition. We consider as a TC model SU(3) gauge theory with 2-\-Nf 
flavors of techniquarks where the mass of two-flavors are fixed to a small value and 
the other TV/ fiavors have arbitrary masses. Nf is also taken to be arbitrary. Based 
on an analysis of linear sigma modef^, we infer that the nature of the chiral phase 
transition at finite temperature changes from crossover to first order at a critical 
mass as the mass of Nf fiavors are decreased from infinity while keeping the two- 
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flavors' mass constant. It should be noted that, if the first order transition is strong 
enough (or equivalently the masses of 7V/-flavors are hght enough), the electroweak 
baryogenesis (EWBG) would become viable. EWBG is, in general, restrictive and 
attractive in that no tunable parameter exists and its success is solely determined 
by the dynamics intrinsic to gauge theory. Then, it is interesting to ask what is 
the upper bound on the masses of 7V/-fiavors which allows the first order phase 
transition and successful EWBG. EWBG in the standard model (SM) was studied 
on the lattice and turned out to fail since mn ^125 GeV is too heavy to induce the 
first order transition.'^ As for TC models, the studies in this direction were carried 
out in the context of the effective theory and obtained promising results.i^ As a first 
step toward the rigorous test of this possibility we study the thermal nature of a 
many-flavor TC model by lattice numerical simulations. To be precise, we aim at 
putting the upper bound on the mass of A^j-flavors of fermions by requiring the 
occurrence of first order phase transition. This upper bound can then be translated 
into that on the technipion mass, which can be directly compared to the results of 
LHC. 

Here let us describe why we consider "2+ A^/" -flavors. In many-flavor TC models 
without the S'C/l(2) x Uy{1) interactions, two flavors of them have to be exactly 
massless and the resulting three massless Nambu-Goldstone bosons (NGBs) are 
absorbed into the longitudinal mode of the weak gauge bosons when one turns on 
the SUl{'2) X C/f(1) interactions. On the other hand, the mass of other Nf flavors 
must be larger than an appropriate lower bound otherwise SxSB produces too many 
(light pseudo) NGBs, none of which is observed yet. Furthermore, in the presence 
of too many massless and almost massless NGBs, S'-parametep^El becomes large or 
even diverges. Thus, A^/-flavors have to have explicit breakings of an appropriate 
size. In this work, we simply let them have a mass. 

Our model based on SU(3) gauge theory is essentially the same as many flavor 
QCD except for their dynamical scales; ~'l TcV for TC and ~^1 GeV for QCD and 
thus we can simply apply numerical techniques developed in lattice QCD to the 
study of TC. As discussed below, the critical mass increases with Nf. Hence, from 
the viewpoint of lattice numerical simulation, the boundary of the flrst order region 
can be reached more easily for large Nf. 

Another purpose of this study is to understand the real QCD with 2-1-1 fla- 
vors. At the physical masses and zero density, the chiral transition is known to be 
crossover, and is expected to become flrst order at a critical density. Toward the 
determination of the critical density, it is important to flnd the critical surface in 
the parameter space spanned by masses and chemical potential! ^^ * ^^ ' However, recent 
lattice QCD studies suggest that the critical surface at zero density is located in the 
very light quark mass region and it makes the determination extremely difficult .1^^ 
Fortunately, some of properties arc independent of Nf. The study of 2+A/-fiavor 
QCD is expected to provide important information for 2-|-l-flavor QCD. 

We flrst describe the method to identify the nature of the phase transition and 
then present the critical mass separating the first order and crossover regions in 
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2+A^/-flavor QCD. The work reported here has been aheady pubhshed in Ref. [T31 

2. Method 

We examine the effective potential defined by the probability distribution function 
of the gauge action to identify the nature of the phase transition. The first order 
transition is concluded by the existence of the two peaks in the distribution func- 
tioniSESl We define the distribution function for 2+A^/-flavor QCD with the quark 
masses nif {f ~ 1,- ■ ■ , 2 + Nf) by 

w{P; 13, ruf) = / VUVi^V^j d{P - P) e^-^'-^" 

„ Nt+2 

= VU 5{P - P) e^^^-'"^ Yl (detMimf)), (1) 



where Sg and Sq are the gauge and quark actions, respectively, and M is the quark 
matrix. iVgito = N^ x Nt is the number of sites. (3 = 6/gl is the inverse lattice bare 
coupling, and P = —Sg/{6NsitcP)- The effective potential is then defined by 

V,s{P;(3,mf) = ~lnw{P;(3,mf). (2) 

We consider QCD with two degenerate light quarks of the mass mi and Nf 
quarks of to/j. For later convenience, the potential is separated into two parts; one 
is the contribution from two-flavor QCD Vo{P; (3) and the other is the rest, 

V,s{P; /3, mh) = Vo{P; (3o) - lni?(P; /3, m,,; /3o), (3) 

with 

n tet mS / ' ^^) 

where (• • • )p:fixed = {^{P—P) ■ ■ ■ )lial {^{P~P))po ^'^'^ (' ' ' ) Po denotes the ensemble 
average over two-flavor configurations generated at /3o and mi . Since the m; depen- 
dence is not discussed in the following, it is omitted from the arguments. /3o is the 
simulation point, which may differ from /?. By performing simulations at various 
/3o, one can obtain the potential in a wide range of P. 

Restricting the calculation to the heavy quark region, the determinant for Nf 
flavors in eq. ^ is approximated at the leading order as 

' det M{Kh) 



In 



detM(O) 
for the standard Wilson quark action and 



288iV,itc4^ + l2Nl{2Khf'n +• • • (5) 



In 



det Af (mh) 



detM(oo) 



36^..^^ 6^^^... (,) 



(2mh)4 (2mh)^* 

for the four-flavor standard staggered quark with TTih- Hh in eq. ([5]) is the hopping 
parameter being proportional to l/rrih, and fl is the real part of the Polyakov loop. 
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Fig. 1. Left; lnR(P;h,0) as functions of the gauge action. Right: The curvature of lni?(P;/i,0) 
for h = 0.01 - 0.07. The circle and square symbols are d?Vo/dP^(P). 



For improved gauge actions sueli as Sg = 6iVsitc/3[co(plaquette) + Ci (rectangle)], 
additional ci x 0{k*) terms must be contained in eqs. ([5]) and ©, where ci is the 
improvement coefficient and cq = 1 — 8ci. However, since the improvement term 
does not affect the physics, we will cancel these terms by a shift of ci. 

Beyond the critical (3 corresponding to the endpoint of a first order transition, 
VeS takes a double-well shape as a function of P, and equivalently the curvature of 
the potential (PVcs / d? P will take a negative value. Observing this behavior usually 
requires a fine-tuning of /3. However, cPVcs/dP^ is independent of /3 and (PVcs/dP^ 
over the wide range of P can be easily obtained by combining data obtained at 
different (3. Thus the fine-tuning of (3 is not necessary in this case-^i^ We focus on 
the curvature of the effective potential to identify the nature of the phase transition. 

Denoting h = 2iV/(2Kft)^' for TV/ degenerate Wilson quarks, or /i = Nf/{A x 
(2mh)^') for the staggered quarks, we obtain 



In i?(P; /3, Kh] /3o) = In i?(P; h) + (plaquette term) 



o{4'^') 



(7) 
(8) 

P:fixcd,^0 

Notice that R{P;h) is independent of Pq. The plaquette term does not contribute 
to d'^Ves/dP'^ and can be absorbed by shifting (3 —i' (3* = /3 + ASNjk'I for Wilson 
quarks. iV/-flavors do not have to be degenerate. The non-degenerate case is realized 



by redefining h = 2J2fiii2Kf)^' or h = (1/4) X;/=i(2to/)-^'. In the following 
we discuss the mass dependence of R through the parameter h. 



3. Numerical results 

We explicitly demonstrate the above method. We use the two-flavor QCD configu- 
rations generated with p4-improved staggered quark and Symanzik-improved gauge 
actions in Ref. [TBI The lattice size A'sito is 16'^ x 4, and the data are obtained at 
sixteen values of /3 from /3 = 3.52 to 4.00 keeping the bare quark mass to ma = 0.1. 
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The number of trajectories is 10,000 - 40,000, depending on /?. The correspond- 
ing temperature normalized by the pseudo-critical temperature is in the range of 
T/Tc ~ 0.76 to 1.98, and the pseudo-critical point is about /3 = 3.65, where the n-p 
ratio is mps/rriY « 0.7. Further details on the simulation parameters are given in 
Ref. [TBI The same data set is used to study the phase structure of two-flavor QCD 
at finite density in Ref. [M] 

We first calculate the potential in two-flavor QCD, Vo{P;l3), the flrst term 
in eq. ([3]). Because the finite temperature transition is crossover for two-flavor 
QCD at a finite quark mass, the distribution function is always Gaussian type. We 
thus evaluate the curvature of Vq using an identity for the Gaussian distribution, 
(PVo/dP^ — 6Nsitc/xp, where xp is the gauge action susceptibility, 

XP = 6N,itc{iP-{P)f). (9) 

The slope of Vq in the heavy quark limit can be also measured using an equation 
derived from eqs. ([3]) and (|4]). When one performs a simulation at /3o, the slope 
is zero at the minimum of Vo{P; /3q), and the minimum is realized at P sa {P)^^. 
Hence, we obtaiiP^ 

'-^^^%^ = -m~MN.... (10) 

The result of d?VQ/dP^ is plotted in the right panel of Fig.[TJ The circles with dashed 
lines are calculated by xp ■ The squares are computed by the numerical differential 
of dVo/dP obtained at the minimum of Vq. dVo/dP are the squares in Fig. [2] It is 
seen that two different methods provide the consistent results. 

In the calculation of R{P; h), we use the delta function approximated by 5{x) sa 
1/(Ay^) exp[— (x/A)^], where A = 0.0025 is adopted consulting the resolution and 
the statistical error. Because R{P\ h) is independent of /3, the data obtained at 
various (3 are gathered as is done in Ref. [M] The results for \nR{P; h) are shown 
by solid curves in the left panel of Fig. [T]for /i = 0.01 - 0.07. A rapid increase is 
observed around P ^ 0.82. It is also important to note that the gradient becomes 
larger with h. 

The second derivative d^ In R/dP^ is calculated by fitting \nR to a quadratic 
function of P with a range of P±0.015 and repeating with various P. The results are 
plotted in Fig. [T] (right), where d'^Vo/dP^ is also shown as the circles or the squares 
with dashed lines. This figure shows that d^ {In R) / dP^ becomes larger with h, and 
the maximum around P ~ 0.81 exceeds d^Vo/dP^ for h > 0.06. This indicates that 
the curvature of the effective potential, d'^Ves/dP^ = d'^Vo/dP'^ - d"^ {In R) / dP'^ , 
vanishes at h ^ 0.06 and for large h there exists a region of P where the curvature 
is negative. We estimated the critical value he at which the minimum of d^Vcs/dP^ 
vanishes and obtained he = 0.0614(69). 

To see the appearance of the first order transition in a different way, we show 
dVcs/dP at finite h for /3* = 3.65 in Fig.[2j The shape of the dVcs/dP is independent 
of /3 because d^Vcs/dP^ is /3- independent. dVcs/dP monotonically increases when 
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Fig. 2. The slope of VcsiP; /3, h, 0) normalized at (^, h) = (3.65, 0) for h = 0.0- 0.1. The squares 
are dVo/dP. 



h becomes small, indicating that the transition is crossover. However, the shape 
of dVcs/dP turns into an S-shape at h ^ 0.06, corresponding to the double-well 
potential. 

We defined the parameter h = 2Nf x {2Kh)^* for the Wilson quark. Then, the 
critical Hhc corresponding he decreases as Khc = [^c/(2iV/)]^/^'/2 with Nj, and 
the truncation error from the higher order terms in k^ becomes smaller as Nj 
increases. The application range of the hopping parameter expansion was examined 
in quenched QCD simulations with Nt ~ 4, by explicitly measuring the size of the 
next-to-leading order (NLO) terms of the expansion.'iS Then the NLO contribution 
turned out to be comparable to that in the leading order at Kh ^ 0.18. Hence, 
this method may be applicable up to around Kh ~ 0.1. For instance, in the case of 



Nf = 10 with Nt = 4, Khc is 0.118. 



4. Conclusion and outlook 

We proposed an agile method to study the thermal nature of many-flavor QCD with 
EWBG in TC in mind, and applied it to the 2+iV/-flavor QCD. Fixing the mass of 
two light quarks, we determined the critical mass of the remaining A^^^-flavors, which 
separates the first order and crossover regions. The critical mass is found to become 
larger with Nf. Further studies using this method are given in Ref. [131 including 
the investigations at finite density. We find that the critical mass increases with /x 
in the 2+Nf-Ravor QCD. 

The next step for the estimation of the baryon number asymmetry in TC scenario 
is to quantify the strength of the first order phase transition. Another interesting 
application of our method is to study universal scaling behavior near the tricritical 
point. If the chiral phase transition in the two flavor massless limit is of second order, 
the boundary of the first order transition region mf (yrih) is expected to behave as 
nil ~ I'^h^'' ~"^h|^ in the vicinity of the tricritical point, (mi, TTih, fj.) = (0, m^^'^O), 
from the mean field analysis. This power behavior is universal for any Nf. The den- 
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sity dependence is important as well, which is expected to be mf ^ |/i|^P^ Starting 
from large Nf, the systematic study of properties of real QCD phase transition is 
possible. 
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